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ABSTRACT 

We revise the physics of primary electron acceleration in the “slot gap” (SG) 
above the pulsar polar caps (PCs), a regime originally proposed by Arons and 
Scharlemann (1979) in their electrodynamic model of pulsar PCs. We employ the 
standard definition of the SG as a pair-free space between the last open field lines 
and the boundary of the pair plasma column which is expected to develop above 
the bulk of the PC. The rationale for our revision is that the proper treatment of 
primary acceleration within the pulsar SGs should take into account the effect of the 
narrow geometry of the gap on the electrodynamics within the gap and also to include 
the effect of inertial frame dragging on the particle acceleration. We show that the 
accelerating electric field within the gap, being significantly boosted by the effect of 
frame dragging, becomes reduced because of the gap geometry by a factor proportional 
to the square of the SG width. The combination of the effects of frame dragging and 
geometrical screening in the gap region naturally gives rise to a regime of extended 
acceleration, that is not limited to “favorably curved” field lines as in earlier models, 
and the possibility of multiple-pair production by curvature photons at very high 
altitudes, up to several stellar radii. We present our estimates of the characteristic SG 
thickness across the PC, energetics of primaries accelerated within the gap, high-energy 
bolometric luminosities emitted from the high altitudes in the gaps, and maximum 
heating luminosities produced by positrons returning from the elevated pair fronts. 
The estimated theoretical high-energy luminosities are in good agreement with the 
corresponding empirical relationships for 7 -ray pulsars. We illustrate the results of our 
modeling of the pair cascades and 7 -ray emission from the high altitudes in the SG 
for the Crab pulsar. The combination of the frame-dragging field and high-altitude 
SG emission enables both acceleration at the smaller inclination angles and a larger 
emission beam, both necessary to produce widely-spaced double-peaked profiles. 
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1. INTRODUCTION 

There has been a steadily growing number of rotation-powered pulsars with emission detected 
at high energies, following the success of the Compton Gamma-Ray Observatory (Thompson 
2001) and detections by ROSAT, ASCA, RXTE and most recently by Chandra and XMM (Becker 
Sc Aschenbach 2002). With the number of pulsars having detected emission above 1 keV now 
approaching several dozen, high-energy radiation seems to be a common feature. Although the 
theory of acceleration and high-energy emission in pulsars has been studied for over twenty-five 
years, the origin of the pulsed non-thermal emission is a question that remains unsettled. The 
observations to date have not been able to clearly distinguish between an emission site at the 
magnetic poles (Daugherty Sc Harding 1996, hereafter DH96) and emission from the outer 
magnetospheric gaps (Cheng, Ho Sc Ruderman 1986, Hirotani Sc Shibata 2001). Complicating the 
problem is the fact that both models have some unresolved difficulties. In the case of polar cap 
(PC) models, while the energetics and pair-cascade spectrum have had success in reproducing the 
observations, the predicted beam size of radiation emitted near the neutron star (NS) surface is too 
small to produce the wide pulse profiles that are observed. However, Arons (1983, hereafter A83) 
first noted the possibility of a high-altitude acceleration region or “slot gap” (SG) near the PC 
rim, based on the finding of Arons & Scharlemann (1979, hereafter AS79) that the pair formation 
front (PFF), above which the accelerating field is screened, occurs at increasingly higher altitude 
as the magnetic colatitude approaches the last open field line where the electric field vanishes. The 
SG is a characteristic feature that is unavoidable in any PC space-charge limited flow model. In 
the last few years, a self-consistent model of PC acceleration has been developed, which includes 
effects of general relativity (Harding Sc Muslimov 1998, hereafter HM98) and screening by pairs 
(Harding & Muslimov 2001, hereafter HM01; Harding Sc Muslimov 2002, hereafter HM02). HM98 
concentrated on acceleration in the interior regions of the PC, but their calculation confirmed that 
a SG also forms near the PC rim in the presence of inertial frame dragging. Interestingly, since 
the magnitude, and longitudinal and radial dependences of the accelerating electric field used 
in HM98 are different from those in the model of AS79, in the calculation of HM98 the PFF is 
curving up more abruptly to form the inner boundary of a much narrower SG. 

The present paper is motivated by the fact that, despite the attractive possibility of primary 
acceleration within the SG region of a pulsar inner magnetosphere (AS79), the electrodynamics 
and energetics of the SG acceleration remain to a large extent unexplored beyond the original 
model of AS79. What we know about SG acceleration is mostly based on a classical model of 
AS79 and further developments by A83 and Arons (1996, hereafter A96), and it has never been 
investigated in any other electrodynamic model thus far. However, we believe that within the 
general-relativistic model of space-charge limited flow from the pulsar PC (Muslimov Sc Tsygan, 
1992, hereafter MT92) the acceleration properties within the SG will significantly differ from 
those implied by the SG originally proposed by AS79. The SG model of AS79 and A83 produced 
effective acceleration only for large inclination angles and only on “favorably curved” field lines, 
i.e. those field lines curving toward the rotation axis. The frame-dragging acceleration dominates 
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at most pulsar periods at small inclination angles, and operates over the entire PC (i.e. it is 
not restricted to favorably curved field lines). Besides this simple update of the underlying 
electrodynamic model, we incorporate an additional effect of screening of the electric field within 
the PC slot, a narrow space between the last open field lines and lines confining the pair plasma 
column generated at relatively low altitudes above the PC surface. This enhanced screening is 
caused by taking into account a second conducting boundary formed by the pair plasma column 
very near the last open field line. The £j| within the SG is greatly reduced by the proximity 
of two conducting boundaries and is well-known from the electrodynamics in narrow tubes and 
cones with conducting surfaces (Landau & Lifshitz 1984). Ironically, this effect of screening 
by two closely spaced conducting boundaries was overlooked in original calculations of the SG 
acceleration, in which only the conducting boundary at the last open field line was taken into 
account. To estimate the field in the SG, A83 simply applied the solution for the accelerating 
electric field in the absence of the pair plasma, assuming screening only from the last open field 
line boundary. 

The major advantage of acceleration within the SG is the possibility of establishing the 
curvature radiation (CR) PFF and generation of 7— ray emission at high altitudes. The need 
to locate the generation of pair cascades and 7-ray emission at high altitude above the PC was 
demonstrated by Daugherty & Harding (1996, hereafter DH96), and later also by Rudak et al. 
(2002) who numerically modeled the pair formation, high-energy radiation spectra and light 
curves at different altitudes in pulsars. They found that the shape of the cascade spectrum and 
high-energy cutoff, due to magnetic pair creation, for the Vela pulsar best matches the observed 
spectrum for emission at altitudes around 1-2 stellar radii above the surface. Furthermore, the 
wide 7-ray profiles observed can be produced in PC models only if the emission beam opening 
angle is comparable to the inclination angle. High altitude emission along flaring field lines would 
enable larger high energy beams. DH96 approached the problem by trying to reproduce the 
observed spectra and light curves of the Vela pulsar by assuming a simple parameterized model 
of the accelerating electric field. However, the physical origin of that accelerating field as well 
as possible regimes of acceleration and energetics for different pulsar spin periods and magnetic 
fields could not be addressed in that approach. Also the details of such modeling, such as the 
pair-formation altitudes, hardness of the spectrum, and widths of the high-energy pulses, will 
also depend on the regime of acceleration of primaries within the SG. In his study of high-energy 
emission from the SG, A96 came to the conclusion that the SG energetics, as determined by a 
potential which did not include the frame-dragging effect, is not enough to account for the observed 
7-ray luminosities and spectra of the Crab and Vela pulsars. Furthermore, he investigated only 
the CR from the primary particles accelerating in the SG which results in a spectrum that does 
not match the softer observed spectra of 7-ray pulsars. In this paper, we additionally consider the 
emission from the high-altitude pair cascades above the PFF that bounds the inner edge of the 
SG. We derive the general relativistic accelerating electric field in the case of the SG to include 
the effect of screening caused by the slot geometry, and present our calculations of pair cascades 
at high altitudes (at a few stellar radii above the surface), 7- ray spectra, and 7- ray light 
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curves. We provide simple expressions for the 7-ray and heating luminosities, and discuss possible 
constraints on the Crab pulsar viewing angles and obliquities. 

The paper is organized as follows. In §2 we discuss the electrodynamics within the SG regions 
of pulsars. We present the electric potential drop in §2.1, calculate the energetics of primaries in 
§2.2, estimate the characteristic gap thickness in §2.3, and derive 7— ray luminosities per solid 
angle emitted by the SGs as well as heating luminosities in §2.4. In §3 we discuss the main results 
of our numerical cascade simulations. Finally, in §4 . we discuss our main results and draw our 
principal conclusions. 


2. Electrodynamics of Slot Gap 

The SG model for the emission region of a pulsar was originally proposed by AS79, who 
were the first to observe that in a space-charge limited flow the process of pair formation may 
significantly vary across the PC. For example, in their electrodynamic model, the pair formation 
surface tends to curve up and become asymptotically tangent to the field lines near the boundary 
of the polar tube in its half with the so-called “favorably curved” field lines. This occurs simply 
because the acceleration of electrons in this region is not sufficient to produce pairs at low 
altitude, leaving a SG around the column of pair plasma. This model was described in detail 
by A83, and also overviewed by Arons (1981) and A96. The very concept of SG acceleration, 
and its potential impact on pair production, radio and high-energy emission from pulsars, was a 
significant finding. However, the original SG treatment failed to produce a viable model of pulsar 
high-energy emission. One reason is that the calculation of electron acceleration within the SG 
was inconsistent with the acceleration and pair formation in rest of polar flux tube region. In this 
paper we discuss and calculate the effect that we believe makes the regime of SG acceleration very 
feasible and attractive. 

Figure 1 schematically illustrates the geometry of the open field lines at the PC. The extended 
SG region is squeezed between the pair plasma column and closed magnetospheric region filled 
with plasma. The ring-like structure on the top indicates the elevated pair formation site. 


2.1. Electric Potential Drop within the Slot Gap 

The SG region is bounded by the last open field lines and the lines with the magnetic 
colatitude (1 — A£s g) times smaller, where A£sg is the latitudinal gap thickness in units of £ 
(£ = 0/0 0 is dimensionless colatitude of a PC field line, and 0 o ~ (QR/cf( l)) 1 / 2 is the footpoint 
of the last open field line). Let us introduce the field line separating the SG into two halves and 
having the magnetic colatitude $o,SG = #o(l - A£sg/ 2), and denote the innermost half-space of 
a SG as region I and the outermost half-space as region II. Then the magnetic field lines within 
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each of these regions can be described by the colatitudes 

r /( 1)1 r 1 

0sg{t)) = @o,sg yj [l - 2 A ^sg^* = 
r /f 1) 1 1/2 T 1 1 

0 ° y 1 - 2 A ^sg( 1 + ^») , ^* € [ 1 , 0 ] for region I ( 1 ) 

and 

r m)i 1/2 i l 

0sg(t?) = <?0,SG |^ + 2 A £sGf* = 

• /( 1) 1 1/2 l I 

0° y 1 ~ 2 ^sg( 1 -^*) , £*€[0,1] for region II, ( 2 ) 

respectively. 

Here £* is the dimensionless colatitude of field lines bounded by the SG regions I and II and which 
is counted from the field line that separates those regions, r] = r/R is the dimensionless radial 
coordinate, and f(r)) is the correction factor for the dipole component of the magnetic field in a 
Schwarzchild metric (see MT92 for details). 

We shall also introduce the SG half-width in the longitudinal direction scaled to the PC 
surface, £sg = (l/2)r pc A£sG, where r pc = ROq is the PC radius. In a small-angle approximation 
(which is perfect for the SG case) the Goldreich-Julian (GJ) charge density and the charge density 
of primary electrons read 

Poi = (i-^cosx+^SG^^sinxcos^e (3) 

and 

P = ~ 2^)5 J(rj C 1 - «)c°sx+ ! 0 sG(l)#(l)sinxcosd> pc , (4) 

respectively. 

Here f) = 2n/P, P is the pulsar spin period, Bo is the surface value of the magnetic field 
strength at magnetic pole, a is the redshift function, k = ( r g /R)(I/I 0 ) ~ 0.15 h 5 /R 3 e is the 
general relativistic parameter entering the frame-dragging effect, r g is the NS gravitational radius, 
I 45 = // 10 45 g-cm 2 , Ro = R/10 6 cm , I is the NS moment of inertia, Io = MR 2 , M and R are, 
respectively, NS mass and radius, if is a relativistic correction factor of order one, x > s the pulsar 
obliquity, and <f> pc is the magnetic azimuthal angle (see MT92 for details). 

For high enough altitudes within the SG, i.e. for 77 - 1 > <5sg , and in a small-angle 
approximation, the general-relativistic Poisson’s equation (see also the derivation of eq. [ 37 ] in 
MT92), 

V ' (a V$ ) = ~ 4rr (P~Po j)> 


(5) 
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translates into 
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( 6 ) 


where 

^4 = k(1 — 1/r} 3 ), (7) 

6 = (3/2)[H(7/)tf S G(>7) - H(1 )^ S g(1)]. (8) 

By employing expressions (1), (2) and changing variables from 8 to £, in eq. (6), and also assuming 
the approximation A£ SG < 1 we get 


Id 2 d 2 \ 

^sg da + d<P 2 pc ) 


$(r), a, 4>\ pc) = -2 $0 <?0 M COS x + B Sin x cos 4 > pc ) , 


(9) 


where 

^SG = 4 = (^) ( 10 ) 

is the screening factor originating from the gap geometry. Also in eq. (9) $o = [SIR/c)Bq R , 
which is the characteristic value of the maximum potential drop in a vacuum solution (see Deutsch 
1955). 


Assuming that the SG surfaces are equipotential ($[£* = 1] = 0), the solution of (9) for the 
electric potential within regions I and II of a SG can be written as 


$(??,£*, d> P c) = *Mo { ^ SG ^(1 -£*)cosx + 2 B 


1 - 


cosh (y/U^a) 


cosh (y^T) J 


sin x cos (j> } 


pc 


( 11 ) 


Since z/ SG <C 1, we can write that 2[1 - cosh ( £, ) / cosh (y/i^r)] ss i/ SG (l - £*). Hence, eq. (11), 
after substitution the explicit expressions for A and B, translates into 


$(*?,£*, </>pc) = <Mfo G * (l “ “3) cosx+ ~ 6 'o,sg^( 1 )x 

(i-a- 


H(rj) / /(l) . 

mrm~ lsinxcos4>pc 


(12) 


Note that the angular size of the SG limits the maximum lengthscale of the angular part of 
Laplacian and naturally enters the solution of Poisson’s equation (see MT92 for details). The 
factor i/ $G , to be evaluated in Section 2.3, is proportional to the square of the SG thickness A£sg- 
Thus, the additional screening introduced by the conducting boundary at the PFF results in a 
reduction of the potential $(»?) and the associated E\\ in the SG. It is important that from the 
electrodynamic point of view, we treat the SG in exact same way we generally treat the whole 
region of open field lines of the PC. Thus, the physics of screening effect within the SG is akin to 
the physics of screening of the vacuum electric field within the open field line region of pulsar by 
the geometry of the region itself. 
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2.2. Primary Energetics within the SG and Emerging 7 - Ray Flux 

Let us write the general formula for the luminosity of primaries accelerating in a SG, 


- , prim( 7 ?PFF) a( - J 


\p(T,mv)dS(v), 


( 13 ) 


■'SG (’’ppf) 

where the integration is over the surface cut in a sphere at the radial distance r] PFF by the SG 
boundaries, and j? ppf refers to the elevated PFF within the SG region. 


In eq. (13) dS is the element of a spherical surface at radial distance 7. After inserting 
expressions (4) and (12) into eq. (13), we get 


£prim — /sG^sd> 


(14) 


where 


/sg — 


SG 


A£ sg j*(l - «) (l - ^3) cos 2 x+ 
H(V) 




1) 



(15) 


V— *7ppp 


is the efficiency of converting the spin-down power into the power of primary particles within the 
SG, and £ s d = fi 4 BfiR 6 / 6c 3 / 2 (1) is the pulsar spin-down luminosity, where Bo/ f(l) is the surface 
value of the magnetic field strength corrected for the “gravitational redshift” . 

Note that the efficiency / SG is a function of pulsar parameters B and P (via B, P - dependence 
of u sa and A£ sc ), so the luminosity of a primary beam scales as L pr ; m oc L“ d , where a = 1/7 and 
4/7 for ££0.1 B cv and ~ 0.1 £ cr , respectively (see § 2.4 for details). 


2 . 3 . Estimate of Slot Gap Longitudinal Thickness 

In this paragraph we estimate the characteristic scale of a SG over the colatitude ^-coordinate. 
We shall discriminate between two cases where the local magnetic field B < 0.1 B CT and B ~ 0.1 B c r , 
respectively, where B cr = 4.413 • 10 13 G is the quantizing value of the magnetic field strength. 
For B < 0.1 £ cr , the attenuation length of the CR photons is determined by the condition 
eB' sin 0kB ~ 0.2, whereas for B > 0.1 B cr it is determined by the threshold condition £sin 0^ > 2, 
where e is the photon energy in units of m e c 2 , B 1 = B/B cr , and 0 i<b is the angle between 
the photon momentum and tangent to the magnetic field line. The explicit formulae for the 
attenuation length read (Harding et al. 1997) 

0.2 p c 


and 


B'e 
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if B < 0.1 B cr 


s = if B > 0.1 B c 
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where p c ^ (4/3) y/cr/Vt is the radius of curvature. 
The case B <0.1 B cr 


Let us estimate the characteristic thickness of the SG, A£ SG , for the case x ~ 0- Since the 
acceleration within the SG occurs mostly in the saturation regime (i.e. at high altitudes, where 
the accelerating E\\ saturates), we can use the following expression for the dimensionless altitude 
zq (here in units of stellar radius) of pair formation due to CR (cf. eq. [40], for the saturated 
regime, in HM01) as a function of £, to estimate the characteristic scale of a SG over £, 


z 0 = 7 x 10 -2 


p r/4 i 

f_OA £ 

B i2 1%' - e? 1 '’ 


(18) 


where Po.i = P/0.1 s, Bu — Po/10 12 G. In our derivation of formula (18) we employed the same 
procedure of minimization of the sum of acceleration length and photon mean-free path, as we did 
in HM98, HM01 and HM02, with the explicit £ dependence of radius of curvature and accelerating 
electric field (for \ ~ 0) included. The PFF is typically relatively flat over the bulk of the PC 
(for 0.2,$£;$ 0.8), and begins turning up as £ approaches the rim of the open field line domain and 
also the magnetic axis. Near the rim the PFF tends to establish at a few stellar radii, so that the 
colatitude £ SG of the beginning of the SG may be estimated via the condition 


dz o 


dt 


'^SG 


A, 


(19) 


where A ~ 0.1 — 0.5 is a free parameter that can be constrained for the known 7-ray pulsars 
by comparison their observed and predicted 7— ray luminosities. The condition (19) means that 
the regime of the SG acceleration sets in when the characteristic variation of the PFF altitude 
with ^-coordinate rescaled by the whole interval of PC £ values (~ 1) becomes comparable to A 
times the stellar radius. Note that formula (19) does not take into account the effect of screening 
near the SG and is based therefore on a somewhat underestimated value of zo which results in 
the underestimation of the SG width. For this reason the most likely value of the parameter A is 
expected to be in the lower range, say A ~ 0.1 — 0.3. Earlier (see HM98 and HM01) we calculated 
the shape of the PFF produced at low altitudes by CR photons. In general, those calculations 
are valid only for the values of £ where the PFF is relatively flat, because they did not take into 
account the effect of electric field screening under consideration in the SG region (which will 
increase the PFF heights there). However, formally, we can extend our calculations for the values 
of £ close to 1, so that we can see how the PFF turns up as we approach the rim. Figures 2a and 
2b show the height of the PFF produced by CR photons at low altitudes as a function of £ near 
the rim region. We present calculations for the values of pulsar spin period P = 0.1 s (Figure 2a) 
and 0.2 s (Figure 2b) and for the discrete values of B, ranging from 0.01 B CT through 0.2 B cr . 
Even though our calculations shown in Figures 2a and 2b are not consistent in the region where 
the PFF is curving up, they may provide a rough idea about the characteristic width of the SG 
region. Figures 2a, b show that the typical values of £ at which the heights of PFF begin to rapidly 
grow roughly agree with our estimates of the SG width (see eqs. [20] and [22] below). 


1 


Assuming that A£ so < 0.3, from eqs (18) and (19) we find that 

A£ sg « 0.2 Po.i(AB 12 )- 4 / 7 / 4 "5 3/7 - 


(20) 


The case B 0.1 B cr 


First of all, note that this case was not treated by HM01, but in fact may be 
high-energy pulsars. Using the threshold condition for the pair attenuation length 
we can write 


z 0 = 7 x 10 -2 


p 7 / 4 i 

MLl £ 


important for 
(see eq. [17]), 

( 21 ) 


so that for the condition (19) we arrive at the estimate, 


A£ sg » 0.2 Po.i(AB} 2 /2 )- 4/7 / 4 " 5 3/7 . 


( 22 ) 


Note that expressions (20) and (22) match at B \ 2 ~ 1 and that A£ SG tends to increase with 
decreasing A. 

By substituting expressions (20) and (22) into formula (12), we can estimate the characteristic 
electric potential drop as 

1/7 

T Volts, (23) 


A$ « 1.5 x 10 12 B{ 2 


where <5 = —1/7 and 3/7 for the low-B and high-B case, respectively, and 


T = 


( 

4 1 \ 

«0.15 \ 

1-3 

V 

TJ 


cos x+ 


3.8 x 10 


-3 


Rl 


1/2 


Po. 1 I 45 , 


/(l) 


H{Tl) \l T] m ~ sinxcos<?i> p c 


(i-a. (24) 


The above formula for A$ illustrates that the accelerating potential (the Lorentz factor of primary 
electrons) is only weakly dependent on the pulsar magnetic field, and the dependence on pulsar 
spin period enters only the second term in expression for 7 which is proportional to sin The 
strong dependence of the parallel electric field on magnetic field and spin period have been mostly 
cancelled by the screening parameter A£ so which depends on the pair production mean-free path. 
The main factors determining A$ are thus, the radius and moment of inertia of a NS, and pulsar 
obliquity angle x (note that Ko.is = hb/R%)- 

By substituting expressions (20) and (22) into formula (23), we can now estimate the 
characteristic Lorentz factor of primary electrons, 7 = eA $/m e c 2 , 
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where 7] acc is the value of rj corresponding to the characteristic acceleration altitude, and 
A 0 .i = A/0.1. For the parameters of a Crab-like pulsar, B \2 = 5, Z45 — 4, and assuming r/ acc = 2, 
and S == —1/7 (low-B case), we get 

7 ~ 3 x 10 7 Ao.f 7 . (26) 

In the present analysis the width of the SG in units of dimensionless coordinate £ is constant 
as a function of altitude. In reality, it may slightly vary with altitude. For example, it may grow 
with altitude because the plasma column boundary may be getting less sharp with altitude, so 
that the efficiency of geometrical screening may slightly degrade (meaning the increase in Ai/ sg ) 
with altitude. Loosely speaking, the effect of fuzziness of the plasma column boundary at high 
altitudes can be attributed to the systematic increase in the deviation of the transverse component 
of the photon mean-free path from the field line toward the magnetic axis. 


2.4. The Gamma- Ray Luminosity per Unit Solid Angle Emitted by the SGs 

Here we calculate the specific high-energy luminosity, L 7 (ft 7 ), radiated by particles accelerated 
within the SGs and emitted into the solid angle f) 7 . 

First, let us calculate the solid angle within which the photons produced by SGs are emitted, 

fi 7 = 2ir / sin Odd, (27) 

-0SG 

min ,7 

where 

(28) 

CL,, ~ "> d l~ 'lr (29) 

Thus, from eq. (27) we can get 

« ^ w 0l V j^A£ SG ster, (30) 

Here r] i 7 is the radial distance of the high-energy emission site in the SG. Using expressions (10), 
(13), (30), and the estimated values of A£ SG (see eqs [20] and [22]), we get 

The case of B ^ 0.1 B cr 

L 7 (0 7 ) = ^i££lfn = 3 x 10 34 e 7 L% 5 P^ 7 R l 6 7/7 A( % ) erg • s" 1 • ster' 1 . (31) 


The case of B ~ 0.1 B cr 
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= 3 X 10 34 ^ Ll^Plf Rl l7 k{ % ) erg • s" 1 • ster" 1 , (32) 


where 


A(t?) 


= 8/7/ 45 6/? {k(! - k) (l - cos 2 x+ 




»(») LM - , 


tf(i)V 7fo) 


n 

sin 2 x j , 


(33) 


and L s d,35 = Xsd/10 35 erg - s'" 1 . In the above expressions the parameter e 1 (0 < e 1 < 1) is the 
efficiency of transforming the primary energetics into the energetics of high-energy emission. 
Note that in the above expressions L 7 (fi 7 ) cx L 3 ^ 7 (low-B case) and oc (high-B case) simply 
because L pr \ m oc L\£ and fi 7 oc L S( (low-B case), and L pr \ m oc and fi 7 a L s ^ 7 (high-B 
case), respectively. 


The above expressions for L 7 (f2 7 ) are equivalent to the observed quantity <J> 7 d 2 , where <2> 7 is 
the high-energy bolometric flux observed at the Earth, and d is the distance to the pulsar. Thus, 
the only free parameters (which may be determined by pair cascade simulations) are the efficiency 
e 1 and the parameter A determining the SG width. All others are intrinsic properties of the NS. 
Figure 3 shows the pulsar empirical (solid circles with error bars) and theoretical (upside-down 
triangles) values of <2 > 7 d 2 as a function of spin-down luminosity, L s( \, The theoretical values are 
calculated for the parameters c 7 = 0.3 and A = 0.1 (see eq. [31]). Note that parameter e 1 can 
range from 0.2 to 0.5 in cascade calculation, and r? 7 = 3. In Figure 3 the dashed line represents 
the absolute upper limit, where the spin-down luminosity is radiated into the unit solid angle, i.e. 
where 4> 7 g? 2 = L s a/1 ster. One can see that there is a good agreement for most high-energy pulsars 
except several of the pulsars, Geminga and PSR B0656+14, having low T s u, and for J0218+4232, 
which is a millisecond pulsar. These pulsars are near or below the CR death line (see HM02, and 
next Section), and therefore have either very wide SGs or no SGs at all. All other high-energy 
pulsars depicted in Figure 3 are above their CR death lines and are expected to have SGs. 


2.5. Energetics of Returning Positrons in Slot Gap 

Let us estimate the maximum possible heating luminosity produced by returning positrons 
within the SG, provided that the PFF establishes at 77 = r] PFF . 

The general equation for the power of returning positrons can be written as (see e.g. MH97, 
eqs.[76]-[78]; and HM01, eq. [61]) 


; *=7 


I (rjppp ) 


*> G jl$d5, 


(34) 
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where the integration is over the surface cut in a sphere at the radial distance rj PFF by the SG 
boundaries, and x + = p + (^ PFF )/|p G j(^ PFF )| is the charge density of returning positrons, scaled by 
the general- relativistic GJ charge density (see eq. [3]). 

Let us calculate, for the sake of illustration, the fractional density of returning positrons for 
X ^ 0 assuming that « I/ 9 gj( 7 ?pff) “ p|/2- Using formulae (3), (4) and representation of k as 
n & 0.15 I^/R% (see § 2.1, right after eq. [4]), we get 

x + ^ 2(i — k/t/ 3 ) ~ ^/^pfp) ^ 2^ ^ ^3 ^ or ^pff ~ (35) 

which means that for a nearly aligned pulsar the maximum fraction of positrons returning from 
the elevated pair fronts of SGs is independent of the PFF altitude, in contrast to the case of the 
PFFs establishing at low altitudes (see HM01, HM02) over the bulk core part of the PC. 

By inserting <f> pc ) and p GJ {rf} into eq. (34) and performing the integration over 0, in 

this case (cf. eqs.[28], [29]), 
from 

^min = Ipff fi 1 )/ /( 7 7pff)(1 “ ^£sg) (36) 

to 

C = ^o^ PFF /(l)//(r/ PFF ), (37) 

we get 

L + = i« 2 (A £ SG ) 3 L sd . (38) 

The case of B <0.1 B cr 


T ^ 4 Y 1 0 3 O K 2 \- 12 /7p* 3 / 7 / 5 / 7 R~~ 6 / 7 T l / 7 pro- 

iU K 0.15 A M).l M5 L sd,35 er § 


.-1 


The case B 0.1 Bcr 


(39) 


L + « 4 x 10 30 k 2 . 15 A- 12/7 P 0 9 ( 7 / 4 5 5 /7 j R- 24/7 L s Y 35 erg 


.-1 


(40) 


The contribution from heating by positrons returning from high altitudes of SGs should be 
taken into account in the estimates of the PC surface temperatures (see HM01, HM02) 

It is interesting to note, that the condition of turning-on of the SG acceleration in pulsars 
should be tightly related to the pair death line condition determined by the pair-formation above 
the central part of the PC (see e.g. HM02 for detailed discussion), simply because the very 
existence of the gap is determined by establishing the pair plasma column and the screening of B|| 
by pairs. The analytic pair CR death line (see eq. [12] in HMZ02) reads 


logP > 


T l°gP-*log/$m -Aj CR| (B,/)-14.6 
| logP - 21og/£t - A ; r CR) (P, I) - 15.4 


P<B 

P~P, 


(CR) 


(41) 



- 13 - 


where 


aJ CR) (P, 7) = 1.5(log/ 45 - 6.3 log Re), 

(42) 

and 


a{ j R) (R,I) = 1.5(log/ 45 - 7.3 log 72$) 

(43) 

are corrections for the NS equation of state, and 


P„ (CR) = 0.1 

(44) 


is the value of pulsar spin period separating the regime of acceleration of primaries where the jE7|j 
is increasing with altitude (unsaturated) from the regime where E\\ remains practically constant 
(saturates). In eq. (41) is the minimum efficiency of converting spin-down luminosity into 
the luminosity of the primary beam needed for pair formation (see HM02 and HMZ02 for details). 
According to our calculation performed in HM02 the value of parameter ranges from 0.2 to 
0.5. For the parameters of a canonical NS and for = 0.2, we get 

logP ~ 2.5 log P - 14, (45) 

as the required condition for SGs to form. 

Note also that the relationship between the pulsar death line and the regime of SG acceleration 
may be justified only for the CR-produced pair front. Even though the SGs may still be formed 
in the case of the ICS pair production above the bulk of the PC, the density of pair plasma may 
not be sufficient to guarantee the efficiency of the geometrical screening. In this case the effective 
parameter Ay sa may significantly differ from that calculated in this paper, including the case 
where it approaches 1 (no geometrical screening within the SG). Our previous calculation of the 
ICS PFFs indicates (see e.g. HM01, HM02) that this situation may well occur in some pulsars 
that are slightly below the death line (e.g. in millisecond and relatively old pulsars). In this paper 
we do not intend to discuss these non-ideal SGs. 


3. Pair Cascade Simulation 

Before we discuss our numerical simulation it is worth mentioning that there is a long 
history associated with the “prototype” of a SG model, the pure geometrical model of pulsar 
emission known as the hollow cone model. By analyzing the smooth variation of the direction 
of linear polarization within the radio pulsar pulse window Radhakrishnan & Cooke (1969) and 
Radhakrishnan (1969) realized that these observations suggested that pulsars must be magnetized 
NSs which produced relativistic particles from near the magnetic PCs. Komesaroff (1970) showed 
that if the radio emission is due to CR then the pulsar beam must be a hollow cone. The double 
pulses observed in some pulsars, the pulse width, polarization, as well as the spectrum of pulses 
found a natural explanation within this model. 
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In order to explore the characteristics of high-energy emission produced by particles 
accelerated in the SG, we have carried out Monte Carlo simulations of the pair cascades initiated 
by primary electrons accelerated by the E\\ in the gap. This calculation is very similar to the 
pair cascade simulation described by DH96. In that paper, primary electrons were accelerated 
above a pulsar PC by a simple model E\\ described by starting and terminating altitudes above 
the NS surface, a constant term and a linear term. These parameters were adjusted to produce 
both a sufficient final Lorentz factor of the primary electrons and a wide enough gamma-ray beam 
geometry to explain the EGRET spectrum and pulse profile of the Vela pulsar. DH96 found that 
a E\\ starting altitude of 2 stellar radii was required, but no physical justification was provided 
to support this value. DH96 also artificially enhanced the primary particle flux at the PC rim 
in order to produce the hollow emission cone, narrow pulses and low level of bridge emission 
necessary to model the Vela profile. 


In the present study, we use the SG E|j — — V4> derived from the potential in eq. (12). The 
acceleration rate as a function of distance is 


^7 

ds 


eEu 


m e c * 


54 (^ COSX+ \ 0ort 1/2 sin X cos <£ pc ) (l-£ 


), 


( 46 ) 


where 5 is defined right after eq. (23). In this formula, for the sake of illustration, we pull out 
the general-relativistic correction from the second term to get its simple explicit ^-dependence. 
Even though in most cases the second term in the above equation is less than the first one, in 
our numerical calculations we use the exact formula for E\\ based on expression (12). The first 
term in parentheses is the contribution to acceleration from the inertial frame-dragging effect, 
where n « 0.15 hs/R^ is the general relativistic parameter that depends on the NS compactness. 
The second term is the contribution from the field-line curvature, similar to the accelerating 
field of AS79. Factors of order one have been omitted from this term in order to simplify the 
expression in the potential of eq. (12). Note that the SG acceleration rate given above is only 
mildly dependent on B and P. We assume that primary electrons begin their acceleration at the 
NS surface on a central field line of the SG, £* = 0. We follow the electrons as they gain energy 
and follow their CR photons through the local magnetic field, accumulating the probability for 
creating an electron-positron pair. If the CR photons pair produce, the synchrotron emission of 
the created pair is simulated and the synchrotron radiation photons are followed in the same way 
to determine if they pair produce or escape. Details of the cascade simulation may be found in 
DH96 and Harding et al. (1997). We ignore inverse Compton emission of the primaries and pairs 
and photon splitting of the CR photons. 


To illustrate the features of the SG acceleration and cascades, we use the period (Po.i = 0.33) 
and surface magnetic field (Bn = 5) of the Crab pulsar for the case B < 0.1 B cr in eq. (46). A 
stellar radius of 16 km is assumed. Although the SG is a region devoid of pair plasma near the 
surface, full pair cascades can develop above the PFFs that form at high altitudes along field lines 
near the PC rim. In contrast to the pair cascades that form in the core of the PC and terminate 
within a stellar radius of the NS surface, the pair cascades from particles accelerated in the SG 
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develop over a larger distance and the screening of the electric field by the pairs thus requires 
a much larger distance. The primary electrons reach Lorentz factors of 10 7 in the SG before 
becoming radiation-reaction limited. However, the pair cascades begin at lower Lorentz factors 
but the pair multiplicity builds up more slowly in comparison to pair cascades in the PC core (see 
HM01), so that the screening distance in the SG is relatively large and can be one or two stellar 
radii. We find that the pair cascade starts at an altitude of 0.5 - 1.0 stellar radii and continues up 
to 3-4 radii. Since the bulk of the escaping high-energy photons are produced between 2 and 4 NS 
radii, the opening angle of the high-energy beam is much larger (about 20 degrees for the Crab 
period) than for near-surface PC emission. 

Figure 4 shows plots of the distribution on the sky of escaping cascade radiation above 100 
keV as a function of phase, </>, and viewing angle with respect to the rotation axis, (, for different 
assumed values of the inclination. For relatively small inclination angles that are comparable to 
the high- energy beam opening angle, emission can be seen over a large range of observer angles. 
As the inclination angle decreases, the hollow cone pattern becomes smaller and emission is visible 
over a smaller range of viewing angles. Also at smaller inclination angles, the emission hollow 
cone is complete because particle acceleration takes place at all phase angles around the magnetic 
pole. This is a result of the fact that the frame dragging term in the potential is independent of 
4 > pc . At larger inclination angles, our results show that the emission occurs over only part of the 
hollow cone. This is due to the sin x c °s <£ pc dependence of the second term in the potential due 
to the “favorably curved” field line acceleration that was the only source of accelerated particle 
in Arons’ potential. At large inclination angles this part of the potential dominates as the frame 
dragging term is declining. Because the favorably curved field lines from the two poles accelerate 
particles in opposite hemispheres, it is impossible for a single observer to view both poles of the 
pulsar in this model unless the high-energy beam opening angle is unrealistically large. Instead, 
at large viewing angles an observer would see a single narrow pulse. 

Pulse profiles are formed by slicing across the phase plots in Figure 4 at a constant viewing 
angle. Examples of pulse profiles at several different inclination and viewing angles are shown in 
Figure 5. It is evident that widely spaced double-peaked profiles can result at relatively small 
inclination by slicing through the middle of the hollow cone. Profiles with two pulses separated 
by more than 100 degrees, as are seen from the Crab, Vela and Geminga, can be formed when 
the inclination angle is comparable to the beam opening angle (about 20 degrees in this model for 
the Crab). Broad single pulses, as seen from PSR B1509-58, may be produced when the observer 
cuts along the edge of the cone. As the inclination angle increases, the maximum pulse separation 
decreases and the pulses themselves also become narrower. Bridge emission as seen in Vela and 
Geminga high-energy profiles (and in the Crab profile around 1 MeV) is not produced in this 
calculation since we have modeled only cascades and emission from the SG, which produces a pure 
hollow cone beam. The core cascades will produce additional emission on field lines interior to the 
SG which will fill in the profile between the two pulses and produce bridge emission. Emission 
from the full PC will be included in future modeling. Emission can be seen in this model at all 
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phases because primary electrons will continue to radiate curvature emission at higher altitudes 
(all the way to the light cylinder) above the termination of the pair cascade. This is the “off-beam” 
emission (Harding & Zhang 2001) that will be much softer than the on-beam emission from the 
cone. 

Our cascade simulation includes aberration of the emitted photons as well as delay times 
between photons emitted at different altitudes. Both of these effects become more important 
at the higher altitudes of the SG cascades and will produce asymmetries in the pulse profiles, 
especially for fast rotators. For example, aberration will produce a pair attenuation cutoff at lower 
energies for the pulse at the leading edge of the profile because the photon emitted at the leading 
edge attain large enough angles to the magnetic field to produce a pair in a shorter distance than 
photons emitted at the trailing edge. This effect was studied recently by Dyks & Rudak (2002), 
who found that the observed disappearance of the first pulse in double-peaked profiles of the Crab, 
Vela and Geminga pulsars could be produced if the photons were emitted at altitudes of several 
stellar radii above the surface. Emission from the SG cascades will therefore naturally produce 
this effect. 


4. Discussion and Conclusions 

In this paper we revise the regime of electron acceleration within the SG regions above pulsar 
PCs. We have demonstrated that our treatment of acceleration within the SG allows us to address 
the longstanding issues of 7-ray production at high altitudes in the pulsar inner magnetosphere, 
interpretation of multifrequency light curves of pulsars, and emission energetics. The fundamental 
differences between the SG model presented in this paper and the original SG model of AS79 are 

1. Our model takes into account the effect of electric field screening within the gap caused by 
the narrowness of the space between the gap equipotential surfaces (effect of nearly parallel 
conducting boundaries on the accelerating electric field). 

2. Our model incorporates the effect of general-relativistic frame dragging, which in most 
cases determines the PC acceleration properties. As a consequence of this, our SG model 
is operable for all azimuthal angles around the PC and for all pulsar obliquities, whereas 
the classical Arons-Scharlemann model operates only for the so-called “favorably curved” 
(curved toward the rotation axis) field lines, and for nearly orthogonal rotators. 

3. We consider the high-energy emission from pair cascades at the inner edge of the SG while 
Arons (1983, 1996) considered only the CR from primaries in the SG. 

The first two effects make it possible to initially boost the electron acceleration and then 
prolong it up to very high altitudes (from few to several stellar radii) where they generate 
pair-producing CR photons. As a result of pair formation at high altitudes, the high-energy 
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photons get emitted into a beam that is larger than the solid angle pertaining to the SG 
itself. These effects also enable the functioning of SGs in all pulsars which are above the CR 
pair-formation death lines (see HMZ02 for the details of our most recent death-line calculations). 
The second property significantly enlarges the allowed parameter space for the pulsar obliquity 
and viewing angles, in particular to small inclination and nearly aligned rotators where the 
original SG model failed to produce any acceleration. Within the SG model proposed in this 
paper both double and single high-energy pulses are possible, and in the case of double pulses the 
separation between the peaks may be anywhere from near zero to 180 degrees. The asymmetry 
of double-peaked high-energy pulses (one peak is higher than other), and misalignment of 7 -ray 
light curves with the X-ray light curves may be attributed to the effect of photon aberration 
which becomes important for high-energy photons generated at high altitudes. Some other effects, 
such as e.g. effect of NS equation-of-state and effects of very strong stellar magnetic field will be 
discussed separately in our subsequent publications. Also, we plan to investigate the details of 
formation of high-energy spectra in SGs of pulsars and perform comparison with the observational 
data. 

In this paper we calculated the predicted high-energy luminosity per unit solid angle from 
the SG of pulsars as a function of pulsar spin-down luminosity. Our theoretical calculations agree 
well with the observational values for 7 -ray pulsars. Given the fact that our model has only two 
major uncertainties, such as the uncertainty in efficiency of conversion of primary energetics into 
the energetics of 7 -ray emission and uncertainty of our calculation of SG width, the agreement 
is rather remarkable. Note that the uncertainty of our SG width calculation translates into the 
uncertainty of order of a few for our calculation of high-energy luminosity per unit solid angle 
which is comparable to a typical uncertainty of the observations. Our calculations depicted in 
Figure 3 were performed for the efficiency of conversion of primary energetics into the energetics 
of high-energy emission e 7 = 0.3, which does not seem to be unreasonable. 

The SG cascades will produce a component of electron-positron pairs at high altitudes, in 
addition to and distinct from the pair component produced by the low-altitude core cascades. 
One might speculate that these two pair components could be the origin of the core and conal 
radio emission components that have been identified from radio morphology studies (e.g. Rankin 
1993, Lyne & Manchester 1988). It will be interesting to explore the properties of these two pair 
components, which may have different characteristics. 

We presented our illustrative calculations of 7 -ray light curves and double-peaked pulse 
profiles for the parameters of the Crab pulsar. We believe that our simulated high-energy light 
curves and pulse profiles can be used to constrain the obliquity and viewing angles for 7 -ray 
pulsars. We plan to perform detailed and extensive simulations of this kind in the future. We also 
plan to perform comparative analysis of two-pole and one-pole models for selected pulsars with 
available multifrequency light curves and indirect studies of their viewing and obliquity angles. 
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Fig. 1. — Schematic illustration of polar cap geometry, showing the outer boundary of the open 
field line region (where E\\ = 0) and the curved shape of the pair formation front (PFF) which 
asymptotically approaches the boundary at high altitude. The slot gap exists between the pair 
plasma which results from the pair cascades above the PFF and the outer boundary. A narrow 
beam of high-energy emission originates from the low-altitude cascade on field line interior to the 
slot gap. A broader, hollow-cone beam originates from the high-altitude cascade above the interior 
edge of the slot gap. is the slot gap thickness (see text) and 0 o ,SG is the colatitude at the center 
of the slot gap. 
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Fig. 2. — Height of the pair formation front for curvature radiation photons above the neutron star 
(NS) surface zo, in units of NS radius R, as a funtion of magnetic colatitude in units of polar cap 
half-angle, £ = 6/9 0 , for different values of surface magnetic field strength B, in units of critical 
field B cr and x = 30°. 







Fig. 3. — Observed flux above 1 keV, <i> 7 , times distance squared (from Thompson (2001)) (solid 
circles) and theoretical values of specific high-energy luminosity from the slot gap, e 7 Lprim/^ from 
eq. (31) (upside-down triangles) vs. spin-down luminosity for known 7-ray pulsars. An efficiency 
of e 1 = 0.3 was assumed. Also A = 0.1, 7? 7 = 3 and the stellar parameters R e = 1.6 and / 45 = 4 
were used. 
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Fig. 4. — Angular distribution of photon emission above 100 keV from pair cascades above the slot 
gap (using a linear 10-level gray scale), for different values of inclination angle x* Theta and phi are 
polar and azimuthal angles with respect to the pulsar rotation axis. Cascade simulations assumed 
values, P = 33 ms, B = 5 x 10 12 G, of the Crab pulsar. Extended faint tracks are high-altitude 
curvature radiation from primary electrons above the cascade region and are an artifact of the 
resolution in magnetic phase angle <f> pc around the polar cap. 
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Fig. 5. — Theoretical pulse profiles of photon emission above 100 keV from the slot gap pair cascade 
at various viewing angles, generated by making cuts through the angular distributions shown in 
Figure 4 at constant phase angle theta. The spikes on the leading and trailing edges of the peaks 
in some profiles are artifacts of the resolution in magnetic phase angle <j) pc around the polar cap. 







